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Abstract 

We discuss spontaneous breaking of internal symmetry and its Nambu-Goldstone (NG) modes 
in dissipative systems. We find that there exist two types of NG modes in dissipative systems 
corresponding to type-A and type-B NG modes in Hamiltonian systems. To demonstrate the 
symmetry breaking, we consider a 0{N) scalar model obeying a Fokker-Planck equation. We show 
that the type-A NG modes in the dissipative system are diffusive modes, while they are propagating 
modes in Hamiltonian systems. We point out that this difference is caused by the existence of 
two types of Noether charges, and Q^: are the symmetry generators of Hamiltonian 

systems, which are not conserved in dissipative systems. are the symmetry generators of 
dissipative systems described by a Fokker-Planck equation, which are conserved. We find that the 
NG modes are propagating modes if Q'^ are conserved, while those are diffusive modes if they are 
not conserved. We also consider a SU{2) x U{1) scalar model with a chemical potential to discuss 
the type-B NG modes. We show that the type-B NG modes have the different dispersion relation 
from those in the Hamiltonian systems. 
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I. INTRODUCTION 


Spontaneous symmetry breaking (SSB) is a universal phenomenon and widely observed 
at various scales in nature, e.g., our vacuum where the electroweak and chiral symmetries 
are spontaneously broken, superconductors, ferromagnets, solid crystals, and so on [1-4], 
Those well known examples are in Hamiltonian systems. 

As is the case with the Hamiltonian systems, it is known that the SSB occurs even in 
dissipative systems such as reaction diffusion [5] and active matter [6, 7] systems. The re¬ 
action diffusion system has translational symmetry. The symmetry is spontaneously broken 
by a pattern structure in space [5] such as the Turing pattern, which is considered as the 
most basic pattern formation in biology [8]. In the active hydrodynamics, which describes 
collective motion of biological organisms such as flocks of birds, the rotational symmetry 
is spontaneously broken. Toner and Tu showed that the Nambu-Goldstone (NG) modes of 
the active hydrodynamics in d dimensions are given as d — 2 diffusive shear modes and a 
pair of propagating sound modes [6, 7]. Recently, the propagating sound mode seems to be 
experimentally observed in the situation where the flocks collectively turns [9]. Attanasi et 
al. also phenomenologically discussed the type of the NG mode, propagating or diffusive 
modes, based on the SSB of the rotational symmetry and conservation law. 

In Hamiltonian systems, NG modes are classified into two types: type-A and type-B [10- 

13] . When a global symmetry Q is broken into its subgroup H, the numbers of type-A (AA) 
and type-B (77 b) NG modes are expressed as 

Na = Nbs - rank([zQ", Q^]), Nb = ^rank([zQ“, Q^]), (1) 

where A^bs = dim(^/"H) is the number of broken symmetries, and Q® are the Noether charges 
(generators) of Q. The total number of NG modes is AA + A^b = A^BS~rank([zQ", Q^])/2 [10- 

14] . For the spontaneous breaking of an internal symmetry, both type-A and type-B NG 
modes are propagating modes. Type-A and type-B NG modes have the linear and quadratic 
dispersion in momentum, respectively. 

Gompared to those in Hamiltonian systems, the relation between the SSB, NG modes and 
their dispersion relations in dissipative systems is not fully understood. We cannot naively 
apply the above argument to the dissipative systems because it is based on the conservation 
law. In the dissipative systems, the symmetry does not mean that the generators of the 
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symmetry are the conserved quantities. For example, a Brownian particle obeying the 
Fokker-Planck equation has rotational symmetry; however, the angular momentum is not 
conserved due to the dissipation and noise. Nevertheless, we know that the SSB occurs 
and there appear gapless modes, in the absence of the conservation law. In the above 
example of the active hydrodynamics, the rotational symmetry is spontaneous broken by 
the expectation value of the nonvanishing velocity held, where the angular momentum is 
not conserved. 

To discuss the symmetry breaking in the dissipative systems, it is useful to consider path 
integral formulation of Langevin or Fokker Planck equation. The formalism is called the 
Martin-Siggia-Rose (MSR) formalism [15, 16], which is successful in the analysis of dynamic 
critical phenomena [17-20]. Even in the dissipative system, we can construct the Noether 
charge Q'^, which is the conserved quantity, as the generator of the symmetry in the MSR 
formalism. This Noether charge is a different charge from that in the Hamiltonian system. 
We refer the Noether charge in the Hamiltonian system to which is not conserved in 
dissipative systems. 

In this paper, we show that the SSB in dissipative systems can be discussed by 
instead of To this end, we will consider a model with 0{N) scalar helds 
When 0{N) symmetry is spontaneous broken into 0{N — 1), nonvanishing order parameters 
{[iQ% (pj^ityX)]) appear. In this case, the NG modes corresponding to type-A NG modes 
become diffusive modes. The existence of diffusive modes is the characteristic feature in 
the dissipative systems because the NG modes associated with spontaneous breaking of an 
internal symmetry in the Hamiltonian systems are propagating modes. We also consider a 
SU{2) X U{1) scalar model with a chemical potential. We hnd that a type-B NG appears in 
the dissipative system. The mode has the quadratic dispersion as in the Hamiltonian system, 
but the damping rate has different order in momentum. In both examples, the absence of 
the conservation of Q'^ is essential to determine the patterns of dispersion relations. We 
nonperturbatively establish these results using the Ward-Takahashi identity for and 

This paper is organized as follows. In Sec. H, to see symmetry in dissipative systems, 
we consider rotational symmetry of Brownian motion as a simple example. We discuss 
the mathematical similarity between the Fokker-Planck and the Schrodinger equations. We 
also give the MSR formalism for readers who are unfamiliar with. In Sec. HI, we consider 
0{N) and SU{2) x U{1) models to discuss the spontaneous breaking of internal symmetries 
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and their NG modes in dissipative systems. In Sec. IV, we establish the resnlts in Sec. Ill 
in the model independent way nsing the Ward-Takahashi identity. Section V is devoted 
to the snmmary and discnssion. In Appendix A, we discuss the action in the real-time 
formalism, which has two types of Noether charges, Q'^ and Q% The action reduces to 
that in a dissipative system by coupling the system to the environment that violates the 
symmetry. 

II. ROTATIONAL SYMMETRY IN DISSIPATIVE SYSTEMS AND MARTIN- 
SIGGIA-ROSE FORMALISM 

We briefly review that symmetry of a Langevin equation can be discussed as in the 
operator and path integral formalism in quantum mechanics [21-25]. To this end, we consider 
the rotational symmetry of Brownian motion. The Langevin equation for the Brownian 
particle x{t) is given by 

^u{t) =-'yu{t) (2) 

where u{t) = dx{t)/dt is the velocity, 7 the friction coefficient and ^{t) the random noise 
that satisfies the fluctuation-dissipation relation, 

(3) 

Here, (...) represents the average over the noise, and T is the temperature of a heat bath. In 
this section, we choose the mass of the Brownian particle unity without loss of generality. 

If ^ and 7 vanish, the system reduces to the Hamiltonian system, and the angular mo¬ 
mentum Ln = X XU is conserved. This results from the rotational symmetry of the equation 
of motion, x ^ x' = Rx and u ^ u' = Ru with a rotation matrix R. The angular mo¬ 
mentum plays a role of the generator of rotational symmetry, = €ijkXk, where 

eijk is the Levi-Civita tensor, and {..., ...}pb represents the Poisson bracket, {xi,Uj}pB = ^ij- 
In contrast, when ^ and 7 exist, the angular momentum is no longer conserved, dL^/dt = 
—'fx xu + xx^^O due to the friction and noise. However, this does not mean that the 
absence of the rotational symmetry. In fact, Eqs. (2) and (3) are still rotationally symmet¬ 
ric under u ^ u' = Ru, ^ ^ = R^. As we will see in the following, this rotational 

symmetry implies that the existence of another conserved quantity. 
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For this purpose, it is useful to introduce the probability distribution of the velocity, 


P{v,t) = {6^^\u{t) -v)). (4) 

The time evolution of P{v,t) obeys the Fokker-Planck equation, 

dtPiv, t) = P(u, t). (5) 

A point is here that we can regard the Fokker-Planck equation as a Schrodinger equation [26]. 
If we rewrite 

d . , 

V ^ q, (6) 

which naturally satisfy the commutation relation [qi-,Pj\ = iSij, the Fokker-Planck equation 
is expressed as 


d,\P{t)) = -H\P{t)) 


( 7 ) 


with the Hamiltonian, 


H = 'jTp^ — ijp ■ q. (8) 

Here, \P{t)) is the state vector, and FT is the Fokker-Planck Hamiltonian. The Fokker- 
Planck equation (5) is identihed as the coordinate representation of Eq. (7). The important 
difference from quantum mechanics is that the Hamiltonian (8) is not hermitian, and thus, 
the left and right eigenstates are generally different. 

Obviously, the Hamiltonian H is invariant under x —)■ Rx and p —)■ Rp, so that we hnd 
the “angular momentum” as the Noether charge of the rotational symmetry, 

La = qxp, (9) 

which commutes with the Fokker-Planck Hamiltonian H. We emphasize that La is not the 
actual angular momentum because La = —iv x (d/dv) ^ Lr. Nevertheless, La plays a role 
of the generator of rotational symmetry; for example, the commutation relation between 
LAi and qj gives 


Q'j] Rijkqk- 


( 10 ) 
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We can also consider the rotational symmetry of states. Let us consider the following state 
l' 0 ( 9 o)) as an example: 


{q\'ip{qo)) = Nexp 


^{q-qo^ 


( 11 ) 


where {q\ is the left eigenstate of q and N is the normalization constant. By operating La 
to \'ip{qo)), we obtain 


{qlLAl'ipiqo)) = -iN{q x go) exp 


^{q-qo? 


( 12 ) 


That is, we have 


T^|' 0 (go)) = 0 for the symmetric state, 

7 ^ 0 for the non-symmetric state. 


as in quantum mechanics. Thus, we can map the Langevin equation to that in “quantum me¬ 
chanics” and dehne the “Hamiltonian” and the “Noether charge.” If the stationary solution 
of the Fokker-Planck equation (7) is a non-symmetric state, the symmetry is spontaneous 
broken. 

By using these, we can discuss the symmetry of dissipative systems and its spontaneous 
breaking as in quantum field theory. 

It is also useful to introduce the path-integral representation of the generating functional 
Z[J] called the Martin-Siggia-Rose (MSR) formalism [15, 26], 

Z[J] = 

= (13) 


where J represents the source. The action S[q,p] is expressed as 

iS[q,p\ = J dt[ip ■ dtq - H] 


= dt 


ip ■ {dtq yg) - yTp 


= 2 / 9 P 


0 


q 

-2jT) \p^ 



(14) 


In the last line, we symmetrized p ■ dtq term as {p ■ dtq — q ■ dtp)/2, and introduced the 
inverse of the retarded and advanced Green functions as 


DnU 


= ±dt + 7 . 


(15) 
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The MSR action is useful to discuss the SSB in the dissipative systems. Namely, we can 
use standard techniques for the SSB even in the dissipative systems as we shall see the next 
section. 

III. SPONTANEOUS SYMMETRY BREAKING IN DISSIPATIVE SYSTEMS 

In this section, we consider SSB in dissipative systems. We hrst provide a toy model 
of scalar helds, (j)°'{t,x), with 0{N) symmetry. We discuss the symmetry and its Noether 
charge based on the MSR formalism. We shall hnd that the Noether charge that we call 
Q'^ is obtained from the symmetry of the action in the MSR formalism. corresponds to 
La of the Brownian particle in the previous section. In addition, if we drop the dissipation 
from our model, another Noether charge arises, which corresponds to the actual angular 
momentum of the Brownian particle, Lji. 

Next, we discuss the spontaneous breaking of 0{N) symmetry in the language of the 
quantum held theory. The corresponding nonvanishing order parameters are {[Q% *)])• 

In our model, we will have all {[iQ% Q^a\) — {[^Q% Q^r\) — 0, so that the NG modes belong 
to type-A NG modes. Furthermore, we show that the NG modes of 0{N) model become 
diffusive modes. This behavior is quite different from those in the Hamiltonian system, in 
which the NG modes become propagating modes, such as the phonon in superhuids. That 
is, the diffusive NG mode is characteristic in the dissipative system. In Sec. IIIG, we also 
consider a SU{2) x U{1) model with a chemical potential, which is a simple model for 
type-B NG mode. In this model, we will have a nonvanishing order parameter, {[iQ% Qr]), 
and thus, there appear a type-B NG mode. We shall hnd the type-B mode has a different 
dissipation relation from that in the Hamiltonian system. 

A. 0{N) scalar model and its symmetry 

Let us consider the following Langevin equation as a toy model; 


StiPiit, x) - x), r}pB = 0, 


(16) 



R 


(17) 
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where 0 ^ and tt^ are the scalar fields that belong to the fundamental representation of 


0{N) symmetry, the subscript a runs 1,2,iV, 7 is the dissipation constant, and {. 
represents the Poisson bracket: 

dX \ f dY \ f dY \ f dX 


■} 


PB 


{X,F}pb = j (fx 


\d(l>R{x)j \d7r%{x) 
is the random noise satisfying 




diT^ix) 


(18) 


,x')) = ASabS{t - t')6^^\x - x'), 


(19) 


where A represents the strength of the noise. If we assume that the fluctuation-dissipation 
relation, A = 2'yT. F is the free energy, which has the following form: 


F[(j)R,'JTR] = / d^X 




( 20 ) 


If we drop 7 and our model reduces to the usual Hamiltonian equation whose Hamiltonian 
is given by Eq. (20). Hence, without the dissipation, our model turns out to be the classical 
version of the Goldstone model that is used to discuss the SSB in quantum field theory [2]. 
The corresponding MSR action and Fokker-Planck Hamiltonian read 


iS= d^x 


dt(i)R - 7r“ 


A 


9t7r^ + 77rR+ ( -V" + m" + m^(0^)" )()>r ) - ^(</>a)^ 


( 21 ) 


and 


H — i'k’aKr + ( 7vr^ + ( —V"' + + -u"'(05j)"' )0^ ) + -^(0^)^) 


( 22 ) 


where tt^ and 0(4 are the canonical momentum, which satisfy the commutation relations, 

[0^(f,a;),<(f,a;')] = [0^(f, a;), 7r^(f, a;')] = 5^^\x - x ), (23) 

and others are zero. The stationary solution of this Fokker-Planck equation is the Gibbs 
distribution. 


{‘t>R,'^B.\P) = P(‘t>R,RR) = 2 ' , 


(24) 


where T = A/{2'y) is the temperature, and Z is the normalization constant such that 
/ d(f)Rd7iRP{(j)R, tir) = 1. \P) is the eigenstate of H with the zero eigenvalue. 












The action (21) is invariant under under the following infinitesimal 0{N) transformations, 


TT^ + Ze“[T"]ab7r^, (25) 

rA^rA+ TT^ ^ TT^+ te^nbA, ( 26 ) 


where e“ is an infinitesimal constant, and is the generator of 0{N) group, which 

satisfies the Lie algebra, [T“,T^] = with the structure constant For example. 


the generator of 0(3) symmetry 

is given as 




^0 0 0 ^ 


^0 0 


O 

1 

. 

o 

= 

0 0 -t 

rp2 _ 

1 

0 0 0 

n~<3 _ 

i 0 0 


o 

O 


y — i 0 Oy 


o 

o 

o 


From this symmetry of the action, we obtain the Noether charges, 

Qa^ f ^A'\T‘‘]ai4>R + ^R'\T‘‘]o.b4>'’A 


(27) 


(28) 


Q\ corresponds to La of the Brownian motion in Sec. II. Q’a satisfies the Lie algebra, 
[Q% ~ If 7 = ^4 = 0, i.e., the Hamiltonian system, the action (21) is invariant 

under another inhnitesimal transformations: 


-)■ 


+ ie“[T"] 




^U, , 

TTa —)■ TT, 


+ ie"[T“] 


2bT^R, 


(29) 


and —)■ 0^, TT^ —)■ TT^. The Noether charges of this symmetry are given as 


CTr^I 


(fx 


^RilrWk 


(30) 


corresponds to the actual angular momentum in Sec. 11. In fact, their Poisson bracket 
satisfies the Lie algebra, iiQ^, Qr}pb = although their commutation relation does 

not, [Qr,Qr\ = 0^. In the MSR formalism, the doubling of the helds, (pR and 0^, occurs 
and for the Hamiltonian system, it causes the doubling of the symmetry. We will further 
discuss the doubling in the realtime formalism in Appendix A. 

We also give the action by integrating and vr^ out: 



^ In quantum field theory, [Qf^, QA may be nonzero. See Appendix A. 
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where d^x = dtd?x and 


= d‘f ± -idt -V^ + m\ 


(32) 


In the next section, we discnss the SSB in this model and the dispersion relation of the NG 
modes. 


B. Spontaneous symmetry breaking in a 0{N) model 

Let ns now discnss the SSB in the 0{N) model with the MSR action (31). When the 
sqnared mass is negative, —= rn? < 0, the symmeteric state (0“) = 0 is disfavored 
becanse the propagator obtained from Eq. (32) contains an nnstable mode, ~ To hnd 
the trne stable state, we look for the stationary solntion of the MSR action (31). For this 
pnrpose, we consider the following potential: 


Kj = m^i4>°A4’R + + «'H4'R?i<P°A4'R, 


(33) 


which is obtained from the action with constant helds. The stationary solntions are given 
by 


50“ 

5Kfr 

50 


— = Z = O5 


(34) 

(35) 


and we obtain the two nontrivial solntions: 


(0 r)^ = and 

o 

/i 


= 0 , 


. 2 /i 2 


(0r) = and (t)\ = 


3kl 


(36) 

(37) 


in addition to 0^ = 0^ = 0 corresponding to the nnstable symmetric state. The second 
solntion (37) is also nnstable, whereas the hrst solntion (36) is stable. To see the instability 
of the second solntion, let us suppose 0^ at a = 1 has the nonzero expectation value: 


/ ia \ ^0 ral 

= ■ 


( 38 ) 


10 




where 0o = h/w. The squared mass for 0^ for b = 2,3,...,N of the second solution is 
negative: 


[D],Ti^ = 0.k = O)= , 




[D-T{u: = 0,k = 0)= - 




0o/V3 


= 0 , 


'eS 




(39) 

(40) 


I0O/V3 

In contrast, for the first solution, we set (0^) = (j)oS°‘^, and obtain the positive squared mass 




= 2fi , and 


= 0 . 


(41) 


<l>0 


4,0 

Therefore, the stable solution is the hrst solution (36) and the SSB occurs. 

As in quantum held theory, the symmetry breaking is characterized by the nonvanishing 
expectation value of Noether charge and a local operator. In our case, it is 


{liQA,4‘l])=ilT“U4‘o- (42) 

There are (A^ — l)’s independent nonvanishing components. Thus, the 0{N) symmetry is 
spontaneously broken into the 0{N — 1) symmetry. The corresponding for b = 2,3,N 
are the Nambu-Goldstone helds. For example, in the case of 0(3), Q\, Q\ and are 

broken Noether charges and NG helds, respectively: 

{[iQ^fp^]) = = 00 7 ^ 0. (43) 


Furthermore, are also spontaneously broken if we take the limit 7 —)• 0 and A —)■ 0: 

([«,0A])=^[nal0O. (44) 

Since ([*Qr,<5^]) = {['>'Q%Qa]) — 0) modes belong to type-A modes. 

Now we discuss the dispersion relation of NG modes. To this end, we consider huctuations 
around the expectation value (38), in which we parametrize the helds as 

(PU^) = (00 + (^r{x), Xr{x)), (p\{x) = {(TaIx), Xa{x)), (45) 

where the subscript b runs 2,3, ...,N. The MSR action in a and x^ turns out to be 
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( 46 ) 












where the inverse propagators ^ and the interaction term i^nt are given as 

± idt - + 2um (47) 

= (48) 

(Ta + 0o(Xr)^ + + XaXr ^20oaR + (xr)^ + • 

(49) 

We can see that x^ does not have the mass term, i.e., x^ is gapless. The dispersion relation 
of the type-A NG modes is determined by 

k) = —— i'^oj + = 0. (50) 

The solutions are given as 

= -\l ± ^^7^ -4fc2 

~ and — *7 H—(51) 

7 7 

where we have expanded u{k) with respect to k up to the second order. Since u{k) has 
no real part, these modes are purely damping modes at small k. One is the diffusive mode 
in which the damping vanishes at fc = 0. The other has a hnite damping even at fc = 0. 
In this model, the number of the diffusive NG modes coincides with the number of broken 
symmetries 77—1. We note that NG modes in Hamiltonian systems become propagating 
mode, such as the spin wave in ferromagnets [13]. In fact, if we set 7 = 0 in Eq. (50), we 
obtain the propagating mode, uj{k) = ±|fc|. As mentioned in the Introduction, the diffusive 
NG mode is characteristic in the dissipative system. 

C. Spontaneous symmetry breaking in a SU{2) x U{1) model 

Here, we discuss the dispersion relation of type-B NG modes in dissipative systems. We 
consider SU{2) x U{1) model with a chemical potential, which is known as a simple model 
for realizing type-B NG modes [27, 28]. Suppose that the MSR action has the form, 

iS = J d^x(^iip\{{-{do + i/if -7 - 75o)(Pr - 2A((p^(Pr)(Pr) 

+ -7 'ydo)ipA - 2\{ip^^ipR))ipA - , (52) 


Gnt — iu 
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where tpi = tpj) is the two component complex scalar helds, fi the chemical potential, A 
the coupling constant, and 7 the friction coefficient. This action is invariant under SU{2) x 
f/(l) transformation, (p, —)■ ^pi + ie^T’^^pi, where T° is the U{1) generator, and T“ (a = 1, 2, 3) 
are the SU{2) generators satisfying the Lie algebra, [T“,T^] = We choose the 

normalization of the generators as trT"T^ = 6°‘^j2. The Noether charges are given as 


Qa = 



(53) 


where tir = {do + ifi)pR and tia = {do + ip — j)pa- We can also dehne 


Q 


a _ 

R — 



(54) 


which are the Noether charges associated with the transformation, pA —t pA + i^aT°‘pR. Let 
us hnd the stationary solution of SS/6pi = 0. Using the same analysis in Sec. IIIB, we find 
a stable stationary solution, pR = (0,u) with v = p/\/^. Since the following expectation 
values. 




(55) 


are nonvanishing, and Q\ corresponds to the type-B NG helds in the nondissipative limit. 
Here V is the volume of the system. To analyze the dispersion relation, we parametrize the 
helds as pr = {x\i + v + 'ip\ + iiP'r), and pA = {x\ + wi) V'i + #i)- Then we hnd the 
inverse propagator for 'ip and y sectors, 

- zyo; + + 2p? 

—2ipoj 

— i^oo + 2ipu 


D/r{uj, k) = 


2ipu 

— ijoj + k^ 


D^R{uJ,k) = 


-2ipu 


-ui'^ — + k^ 


(56) 

(57) 


The dispersion relation is obtained from det = 0 and det = 0. At small k, we hnd 
the dihusive NG mode, u = —i\k\‘^/'y in the ip sector and type-B mode, 

in the y sector. In the limit of 7 —)■ 0, we obtain the dispersion relations ui = ±\k\/y/3 and 
u = ±|fcp/(2/r) in the ip and y sectors, respectively. Therefore, in this model, the type-B 
NG mode is the propagating mode with quadratic dispersion, while the type-A NG mode is 
the dihusive mode. 
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IV. WARD-TAKAHASHI IDENTITY IN DISSIPATIVE SYSTEMS 

In the previous section, we discussed the NG modes and their dispersion relations associ¬ 
ated with the spontaneous breaking of 0{N) and SU (2) x 17(1) models in the classical approx¬ 
imation. In this section, we nonperturbatively establish the result using the Ward-Takahashi 
identity. We consider a system describing a Fokker-Planck equation, dt\P) = —H\P). We 
assume that the Fokker-Planck Hamiltonian does not explicitly depend on time, dtH = 0, 
and the real part of right eigenvalues of H are non-negative and it contains at least one 
zero eigenvalue. In general, the stationary state with the zero eigenvalue may not be the 
thermal state, i.e., a steady state is allowed in this formalism. We assume that the station¬ 
ary state does not break the spacetime symmetries. We consider a continuum symmetry 
group Q with the generator Q°^ as internal symmetry, which commutes with the Fokker- 
Planck Hamiltonian, [i7, Q‘^] = 0. For helds belonging to a linear representation, 0^ and 0^ 
transforms 

[iQ% (I^r] = i[T'"]ab(l)R, and 0^] = i[T“]„b0^. (59) 

We assume that the Poisson bracket is dehned, and we also introduce Qr such that 

( 60 ) 

In general, does not commute with the Hamiltonian, but it does in the Hamiltonian 
system. When we consider the system in which the Poisson bracket cannot be dehned, we 
need not consider in the following argument. 

In the following, we show the following relation from the Ward-Takahashi identity, 

[D-‘]g'(w = 0,fe = 0)([!«, ,#.“]) = 0, (61) 

where DP is the inverse propagator and indices i and j run R and A. The derivation in 
this section is commonly used in quantum held theory [29]. 

To drive Eq. (61), it is useful to move to the path integral representation of the generating 
functional: 

Z[J] = J (-02) 

where we used the vector notation: cf) = (0^,0^) and J = {JAiJpd- assume that 0“ 
contain the order parameters. The action S[4>] and the path integral measure P0 is invariant 
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under the infinitesimal transformation, 0 —where ea is an infinitesimal constant, 
and S'X4>= [iQ'X, 4>] in the operator formalism. Since the generating functional is invariant 
under the reparameterization of the fields 0 —?■ 0' = 0 + €^(520, we have 


Z[J] = J T:)(f)'e^sm+^!d^xj;p' 

= [ + lea [ d^xJ ■ 5^0 ) + 0{e^^ 


= Z[J] (^1 + ze« j d^xJ ■ {5l4>)j^ + 0(6^), 


where and is the expectation value in the presence of the source J: 

1 


(0)j = 


Z J 


From Eq. (63), we obtain the identity 


j d'^xJ ■ {[iQ%(f)])j = 0, 

where we used {S'X<p)j = {[iQ% 0])j- Using the effective action 


r[ 0 ] = —ilnZ[J] — J d^xcl) ■ J■, 


(63) 


(64) 


(65) 


( 66 ) 


we can write the identity as 

/ Jit ^ 

where we used 5T/5cj) = — J. Differentiating Eq. (67) with respect to and taking 

the limit of J —?• 0, we obtain 

j d‘^x[D-^]‘^^{t -t,x - x)(\iQ% 05(t, a;)]) = 0, (68) 

where the inverse of propagator is obtained as t, x —x) = ((5^r/50“(t, x)6(j)’j{ti x')) 

In momentum space, we obtain Eq. (61). This identity represents the eigenvalue equation 
with the zero eigenvalue, whose eigenvectors are ([zQ^,0^]). The number of independent 
eigenvectors is equal to the number of broken Noether charges. 

When is conserved, we obtain a similar result for 


[D-‘]“> = 0,fe = 0)([j<3S,,#.‘]) = 0, 


(69) 
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For 0{N) model with the expectation value (0^) = ^“^00, using useful normalization of 
Q% we can write {[iQ% for a = 2, 3,..., iV. Then, from Eq. (69), we obtain 


[B-‘]S(i^ = 0.fc = 0) = 0, (70) 

This identity gives the constraint to the dispersion relations. 

Now, we expand the inverse of propagators with respect to u and k as 

-ilD-XK = 0, (71) 

-AoXar = S‘"‘C,0R)k^ - trXim^R - + ■ ■ ■ . (72) 

(73) 


where the coefficients C[n,m) and A(^n,m) are generally nonzero without constraint from an¬ 
other symmetry. Here, generally vanishes due to conservation of the proba¬ 

bility. Furthermore, C(o^o) becomes zero from the constraint Eq. (70). From the equa¬ 
tion [D' k) = 0, we obtained the dispersion relation for the dissipative NG modes 

oj ~ —i\k\^. In addition, when, is also conserved, which corresponds to the Hamil¬ 
tonian system, we hnd that ^4(0,0) vanishes from Eq. (69). When the system satishes the 
fluctuation-dissipation relation, H(o,o) is related to 24(o,o) = 2C'i,o7" with the tempera¬ 

ture T, and thus, C'(i^o) vanishes. In this case the NG modes are the propagating modes with 
the dispersion relation u = aR\k\ — m/lfcp [13], where gr and a/ are constants depending 

The symmetry breaking pattern discussed in the 0{N) model for iV > 4 is relatively 
a simple case because the broken Noether charges transform as the vector representation 
under unbroken 0{N — 1) symmetry. In other words, the NG modes belong to the vector 
representation of the 0{N — 1). The unbroken symmetry restricts coupling between NG 
modes and others, and thus, the inverse of the propagators proportional to the Kronecker 
delta^. If the diffusive NG hied has no internal unbroken charge, the analysis will be more 
complicated. In particular, the coupling between NG modes to hydrodynamic modes must 
be taken into account. 

A different type of dispersion relations will be found when {[Qa^Q’r]) is nonzero, which 
corresponds to Type-B modes. For SU{2) x U{1) model, the unbroken symmetry is G(l); 

^ When the unbroken symmetry is the antisymmetric tensor 0(2), there is a possibility to have in the 
inverse of the propagator. 
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we have two second rank invariant tensors in the real representation: the Kronecker delta 
5“^ and the antisymmetric tensor Then, the inverse propagator is expanded as 

-i\D-XR = + ... . (74) 

In this case, the dispersion relation has the form, u = — ic/lfep, where cr and c/ 

are constants. Therefore, the type-B NG mode can propagate. We note that the dispersion 
relation of type-B NG modes in the Hamiltonian system is u = where 

and bj are constants [13]. 

For more general cases, the coefficients are matrices, 

k) = + ■ ■ • . (75) 

If detG[\^Q^ is nonzero, C^ 2 o) is negligible at small k, and we obtain the dispersion relation 
from the eigenvalue of The eigenvalue is generally complex, and thus, we 

have the form u = — idi\k\‘^. We expect the existence of the coefficient is related 

to the nonvanishing {[Q% Q/jj), although we have not given the proof, which is beyond the 
scope of this paper. 

V. SUMMARY AND DISCUSSION 

We discussed spontaneous symmetry breaking (SSB) and the Nambu-Goldstone (NG) 
modes in dissipative systems described by a Langevin or Fokker-Planck equation. For 
this purpose, we employed the 0{N) and SU{2) x U{1) scalar models as toy models. In 
the nondissipative limit, which corresponds to a Hamilton system, there exist the Noether 
charges Q'^ that are the generators of the internal symmetry by means of the Poisson bracket 
and they are conserved. In contrast, in the dissipative system, Q'^ are no longer conserved 
due to dissipation and noise. Instead, there exist other conserved quantities Q% which are 
the Noether charges of the internal symmetry in the Fokker-Planck equation. 

Symmetry breaking is characterized by existence of a nonvanishing order parameter. In 
the 0{N) model, 0{N) symmetry is spontaneously broken into 0{N — 1), and the order 
parameter is {[iQ% ((>%])■ Since {[iQ% Qr]) = 0, the NG modes belong to the type-A modes. 
We found that the NG modes are the diffusive modes, u ~ —ik!^. This is the different 
behavior compared with the Hamiltonian system, where the NG modes are the propagating 
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modes. This difference is caused by whether is conserved: When both Q\ and are 
conserved and their symmetry is spontaneously broken, there appear the propagating NG 
modes. In dissipative systems, when Q'\ is broken, the diffusive NG modes appear. We 
established this result by using the Ward-Takahashi identity for Q\ and Qr symmetries in 
Sec. IV. 

We also discussed type-B NG modes in SU{2) x 17(1) model, where ^ 0. In 

this case, the dispersion relation of NG modes have the form of a; = a|fcp — ih\k\^, while 
they are u = a'\k\^ — h'\k\^ in the Hamiltonian system, where a, b, a', and b' are constant 
parameters. In contrast to type-A NG modes, type-B NG modes are still propagating modes. 

In this paper, we focus only classical systems. Generalization to quantum systems is 
straight forward: We may add higher terms in 0 a such as take into account 

the Bose and Fermi statistics. 

Our approach can apply to the spontaneous breaking of spacetime symmetries, although 
our result in this paper is limited to that of internal symmetry, although the situation will be 
more complicated. Even in the Hamiltonian system, the general counting rule and dispersion 
relation of their NG modes have not been well-known. An interesting example in a dissi¬ 
pative system is discussed in the active hydrodynamics, where energy and momentum are 
not conserved, but equations of motion respect spacetime translational and rotational sym¬ 
metries [6, 7]. In this situation, the velocity fields is the order parameter and it breaks the 
rotational symmetry. For d-spatial dimensions, there appear d — 2 diffusive (shear) modes, 
and one propagating sound mode. This sound mode caused by the mixing between longitu¬ 
dinal NG mode and the hydrodynamic mode associated with the number conservation [6, 7]. 
The mixing of hydrodynamic mode can change the dispersion relation. 

It is interesting to clarify the relation between the broken symmetry, the NG modes 
and their dispersion relations in dissipative systems. Recently, in the Hamiltonian system 
without the Lorentz invariance, those relations have been made clear [10-13], which does not 
cover the dissipative system. From our observations, we propose the following conjecture: 
We suppose that the Fokker-Planck Hamiltonian H commutes with the generator Q°^ of a 
Lie group Q. We also suppose that the Poisson bracket is defined, and exists as the 
generator of Q in the sense of Poisson bracket. In general, does not commute with the 
Fokker-Planck Hamiltonian. When the Q is spontaneously broken into its subgroup 77, the 
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number of type-A (A^a) and type-B (A^b) NG modes will be given as 


(76) 

(77) 


jVa = iVes - rank{|!(3;j, (jJ]), 


where A^bs = dim(^/'H) is the number of broken symmetries. These equations correspond 
to Eq. (1) for the Hamiltonian system. Their dispersion relations will be classified as four 
types: 


Type-A 


Type-B 


uj = cji\k\ — icj\k\‘^, if [Qr, = 0, 
u =-ici\k\'^, if [Q^,iJ] 7 ^ 0, 

00 = Cijlfcp — if [Q% H] = 0, 

00 = Cfilfcp — ic/|fcp, if [Q% H] 7^ 0. 


(78) 

(79) 


Of course, the models discussed in this paper satisfy these relations. We leave the detailed 
analysis and the proof of this conjecture leave to our future work. 
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Appendix A: Symmetry of Hamiltonian and dissipative systems in the real-time 
formalism 


Here, we discuss symmetry of an action in Hamiltonian and dissipative systems from the 
real-time formalism. By coupling the system to the environment explicitly violates the 
symmetry, we shall see that our model corresponds to a dissipative system. 

We consider the following Lagrangian to discuss the symmetry of the Hamiltonian system. 


c = lidtrr - - Imurr - 

The generating functional in the real-time formalism is given as 


(Al) 
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FIG. 1. The contour on the complex time plane in the real-time formalism. The branch Ci runs 
on the real axis from —oo to oo, while C 2 runs backward from 00 to — 00 . runs from —00 to 
—00 — i/3, where j3 is the inverse temperature. 

with the action, 

02] = - ^ J d^xd^x'(j)i{t, -t\x- x) 

d^x(^{{(j)^{t,x)fy - (A3) 

where the subscript i,j run 1,2, and 0^2 represent the fields on the forward and back¬ 
ward branches, respectively (See the caption of Fig. 1.) Here, are the inverse of the 


propagators, which are given as [30] 

Dn{uJ, k) = iP— - ^ - 2 + (b + k), (A4) 

D22{(x, k) = {Dii(u, k))*, (AS) 

Duiuj, k) = n{u)p{u, k), (A6) 

D 2 i{uj, k) = {1 + n{u))p{u, k), (A7) 


where P denotes the principal value, n{u) = l/(e^^ — 1) the Bose distribution function, and 
p{u,k) = 27re(uj)d(uj'^ — k^ — ml) the spectral function with the sign function e{u). The 
important point is that the doubling of the fields occurs: 0“ —>■ 0“ and 02. This doubling 
causes the doubling of the symmetry. 




(A2) 
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Before discussing the symmetry, we change the held variables to 


4>‘r = lir, + 4>i). 4 >‘‘a = 4‘\- (A8) 

The generating functional and the action for 0^ ^ is written in the new variables as 

iS[(j)R, 0 a ] = - ^ J d^xd^x(j)i{t, x)D^^{t -t\x- x)(j)‘^{t\ x') 

- iul j , (AlO) 

where the subscript i and j run R and A, and J^ = and = (J“ + T2)/2- Here, 

D]^\ is given in momentum space as 


0 


k) 


D^\u,k) =1 . ^ V ^ . I 

Wj^\u,k) + n{u)YDj^\u,k) - D/{u,k)) j 

where k) = —{u + ie^ + + ml and D^^{u, k) = k))*. D^\ is inhnitesi- 

mally small in this case because 

'1 


Da\ = ^(2 + i^R ~ 

= 2ea;(l + 2n{u)) ~ e. 


(A12) 


Let us discuss the symmetry of Eq. (AlO). We can easily see that the action is invariant 
under the following transformation: 


« ^ « H- te“lT“UA, + k“[T“UA- 


(A13) 


The Noether charge of this symmetry is given as 


Q'X = / d^x 


Tr‘\i[T%b(l)R + T^Ri[T%b(l)A 


(A14) 


where the tira are dehned as 


"■r = ("■“ + "‘ 2)/2 , JTr = - rS- 


(A16) 


Here, vr“ 2 = 2 are the canonical momentum of 0“ 2 and satisfy the commutation relations 

[<Pl{t,x),7rl{t,y)] = - y), a;), 7r2(t, ?/)] = -i6^’’6^^\x - y), (A16) 
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We note that 02 7rf are the helds on the backward branch C 2 hence the commutation 
relation has the negative sign. We can see that the form of the charge (A14) is equal to that 
of the Langevin equation (28). 

Furthermore, in the limit e —?• 0, Eq. (AlO) is invariant under the transformation, 

« ^ (A17) 


and its Noether charge is given as 


Q 


OL _ 

R ~ 






(A18) 


We note that this charge corresponds to the charge (30) if we drop {1 /4:)7i'\[T°‘]ab((>A term 
in Eq. (A18). Therefore, the action of the closed system (AlO) is invariant under the 
transformations by and Q^. That is, the doubling of the symmetry occurs. In the original 
Lagrangian (Al) is invariant under the transformation 0“ —)■ 0“ + ie“[T"]afe0^. Meanwhile, 
the action of the real-time formalism (AlO) is invariant under the two transformations, 
Eqs. (A14) and (A18). 

Next, we consider the system of 0“ couples with other environment scalar helds with 
mass MK We assume that the interaction between 0“ and has the following form: 




(A19) 


where g is the coupling constant. Then the total action is S'open = >S'[0“]-|-S'[<I)^’“]-|-S'int[0“, 
with 




0 iG' 




4> 


l,a 

R 




^l,a 

A 


(A20) 


where 


—-{u ±ieY + + Mf, — i(^-+ (A21) 

in momentum space. 

By integrating the environment helds out, we obtain the ehective action for 0“: 

- ig'^Gj 


iSeS = -X 


du dk’k 
^ (27r)3 

o f duj dkk 


0 - ig-(JA 

.2/ 




iDjR - ig^Gn - ig^G 


rR 


lUn 


2-k (27r) 




(A22) 
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where Gr^a = Yhi ^R,Ay 



(A23) 


We take continuum limit with respect to the index I of G^^^, and replace 



'OO 



(A24) 


where is a weight function. Hence, Gr turns out to be 



(A25) 


Now, we discuss the symmetry of the action Eq. (A22). The important point is that it is 
not symmetric under Eq. (A17) because Gxih) ^ 0; see Eqs. (A22) and (A25). That is, the 
Q'r symmetry Eq. (A17) is explicitly broken. Conversely, symmetry remains even in the 
dissipative system. The original action Sopen is invariant under the simultaneous rotation 
0“ and whereas the integrated action Eq. (A22) does not have fields. As a result, 
symmetry is broken, while Q\ symmetry remains as the remnant of the original symmetry. 
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